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I. INTRODUCTION 


This report concerns our research conducted on Lake Washington in collaboration with S.P 
Gogineni and R.K. Moore of the University of Kansas. The aim of the research was to investigate the 
relationship between wind, waves and radar backscatter from water surface. To this end, three field 
experiments with periods of 2 to 4 weeks were carried out during summer months in 1988, 1989 and 
1990. For these periods, the University of Washington group provided (j) environmental parameters 
such as wind speed, wind stress and atmospheric stratification through measurements of surface 
fluxes (of momentum, sensible heat and latent heat) and of air and water temperatures; (u) wave 
height spectra including both the dominant waves and the short gravity-capillary waves. Surface 
flux measurements were performed by using our well tested instruments: a K— Gill twin 
propeller-vane anemometer (Ataktiirk and Katsaros, 1989) and a fast response thermocouple 
psychrometer. Wave heights were measured by a resistance wire wave gauge described by Ataktiirk 
(1991) and Ataktiirk and Katsaros (1987). The University of Kansas group was responsible for the 
operation of the microwave radars. 


II. ACCOMPLISHMENTS UNDER THIS GRANT 

The accomplishments and highlights of the experimental results are summarized below. 
i. Surface Fluxes and Sea State: 

Many attempts have been made to parameterize the surface fluxes in terms of routinely measured 
bulk quantities. A major difficulty in these attempts has been characterizations of the poorly known 
surface roughness parameters, z 0 , z, and z^ for wind speed, sensible heat and water vapor, 
respectively (Blanc, 1985, 1987). 

At our field site on Lake Washington we have the proper exposure for studying the surface fluxes of 
momentum, sensible heat and water vapor over pure wind waves at wind speeds of less than 10 m/s 
and a fetch of 7 km. Analyses of our Lake Washington data collected over the past few years show 
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that z, and z q do not vary with wind speed (stress ) or sea state. Within the wind speed range of 3.5 to 8 
m/s encountered in the field work, z, and z () are approximately constant. On the other hand, z„ 
depends on the sea state such that 
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where C p is phase speed of waves at the peak of the frequency spectrum, u* is friction velocity and g 
is acceleration due to gravity. This result based on the atmospheric turbulence and wave 
measurements is consistent with other results obtained from field experiments (Merzi and Graf, 
1985; Geemaert et al., 1987; Donelan, 1990; Smith et al., 1992). 

We also devised a technique to estimate z 0 from rms wave height and the equilibrium range 
parameter of the measured wave spectrum. For the cases of pure wind waves subject to unsteady 
wind speeds on Lake Washington, the technique provided better estimates of z 0 hence, the neutral 
drag coefficient, than the atmospheric turbulent flux measurements. If the response of the waves to 
turning winds is included, the approach may have the potential to determine the momentum flux in 
frontal zones where the only direct method (i.e. eddy correlation technique) of measuring surface 
fluxes is no longer applicable due to lack of stationarity and horizontal homogeneity required by the 
similarity theory. 

The results and relationships given above are analyzed in the Ph.D. dissertation by Ataktiirk ( 1991 ) 
and were presented (Ataktiirk and Katsaros, 1991 ) in Air-Sea Fluxes and Water Waves Symposium, 
XX General Assembly of the IUGG. 11-24 August 1991, Vienna, Austria. 

ii. Wave Breaking Statistics: 

By analyzing video films obtained simultaneously with our wind stress and wave height data, we 
also studied the incidence of breaking waves of three kinds: microscale, spilling and plunging types. 
Relationships were sought between various environmental parameters and the frequency of 
occurrence of each class of breaking. 

Frequency of occurrence of microscale breaking was uniform in all the cases studied. However, 
those of spilling and plunging breakers showed a clear dependence on the wind stress and the 
evolutionary stage of the water waves hence, on the atmospheric turbulence intensity. Regression 
analysis indicated that percentage of crests with spilling and plunging breaking, %Bc + p, could best 
be predicted by 


%B S+ = 32 . 87 ul + 42 . 37 ^ - 2 . 47 
1 t p 

with a correlation coefficient of 0.66 and standard error of 0.93. 
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Interpretation of the above relationship is straightforward. Since breaking results from excess wave 
energy, %B is expected to be larger for higher u* 2 (i.e., wind stress). Also, for a given u*. breaking is 
more frequent in a developing wave field (smaller C p , larger u*/C p ) than in a mature wave field 
(larger C p , smaller u*/C p ). 

The wind input received by the waves is distributed among the wave components through non-linear 
wave-wave interactions. The temporal and spatial scales of this process are closely related to the 
characteristics of the wave field. Therefore, we hypothesize that the absolute width of the wave 
spectrum (or the value of the wavenumber at the peak of the spectrum) may be an important factor in 
wave breaking statistics. This subject requires further study. 

We were also able to relate the variance of wave height in acertain frequency band. 6- 1 0 Hz. with the 
wave breaking events. The visible evidence of breaking as spillers orplungers occurred for variances 
above a certain threshold. Therefore, wave height measurements provide an objective means of 
identifying wave breaking events. 

The above findings of Katsaros and Ataktiirk (1991) were presented in IUTAM Symposium on 
Breaking Waves, 15-19 July 1991 , Sydney, Australia. 

iii. Slicks - Effects on the Wave Spectrum: 

The effects of surface slicks, man-made and natural, on different scales of waves in the spectrum 
were also studied with these data (Katsaros et al., 1989; Ataktiirk, 1991). The results showed that 
surface slicks of either origin caused suppression of small scale waves in a similar manner. Waves 
with length of about 5 cm were found to be affected the most in accordance with the Marangoni 
theory. Influences on the dominant waves as observed occasionally were hypothesized to be due to 
combined effects of surface slicks and enhanced intermittency in atmospheric turbulence during 
conditions of stable atmospheric stratification. 

iv. Radar Backscatter Studies: 

In parallel with our hydrodynamic studies we collaborated with colleagues from the University of 
Kansas, Professors P. Gogineni and R.V. Moore, and their students. Simultaneous measurements of 
atmospheric turbulence, waves and radar backscatter power and range (at X and C bands in 1 988 and 
1989, and at Ka band in 1990) were conducted at our site on Lake Washington. 

Preliminary processing of the 1989 and 1990 radar data by the University of Kansas group has been 
completed. The results include (i) time series of backscattered power, wave height and normalized 
radar cross section, (») power spectral densities of radar backscatter and wave height, and (iii) 
magnitude and phase of the modulation transfer function for various combinations of polarization, 
incidence angle and up/cross wind look angle (Salam et al.. 1991). 

The University of Kansas group also carried out a study of the sea spikes in the radar backscatter 
signal as a function of wind stress (Bush et al., 1 991 ). Their results showed that while some of the sea 
spikes were associated with breaking waves or wedges (steep waves with no spilling ) . most of the sea 
spikes (about 75%) could not be attributed to any visible surface features. 
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On our part, analyses of the atmospheric environmental variables and surface waves corresponding 
to the periods of joint efforts were completed. Other plans for the current funding year include joint 
consideration of the available data sets to relate the radar backscatter signal to environmental 
variables, to compare statistics of breaking waves determined independently from radar backscatter 
and wave height measurements and, to investigate the modulation of short waves bv the underlying 
long waves. 

ID. TECHNICAL REPORTS, PRESENTATIONS 

One of the major accomplishments under grant NAGW-1 322 to the University of Washington from 
the Ocean Processes Branch of National Aeronautics and Space Administration is the Ph.D. 
dissertation of Dr. S. S. Ataktiirk (Ataktiirk, 1991 ). Copies of this thesis were provided to NASA as a 
Technical Report. 

The results from the above study were presented in Air— Sea Fluxes and Water Wtoi rw Symposium. 
XX General Assembly of the IUGG, 11-24 August 1991, Vienna. Austria (Ataktiirk and Katsaros. 
1991 ). Following this symposium, we participated in a workshop to review the state of knowledge of 
dependence of drag coefficient on wave properties. The workshop was concluded with a statement 
of where we are and what needs to be done on this subject (see Appendix A). 

Our findings on the dependence of wave breaking statistics on wind stress and wave development 
were described in a manuscript by Katsaros and Ataktiirk (1991). Copies of this manuscript were 
provided to NASA. (An additional copy is also enclosed with this report as a reference.) 

The results from wave breaking study were presented in the International Union of Theoretical and 
Applied Mechanics, IUTAM, Breaking Waves Symposium, 15-19 July 1991 , Sydney. Australia. 
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Appendix A: 


FROM: Kristina B. Katsaros 

TO: Prof. Yoshiaki Toba, Sendai, Japan 

Prof. Ian F. Jones, Sydney, Australia 


Message: Ideas for the summary statement of Workshop on Surface Fluxes and Water Waves 

Wind/wave research, in spite of dedicated efforts by numerous very talented scientists, is still in early 
stages of development. The physical-dynamic processes involved in wave growth are not 
sufficiently well understood that we can parameterize wave development in terms of a few 
parameters for most of the complex and transient situations where several wave trains exist 
simultaneously or the synoptic situation is not steady. However, progress has been made such that the 
simplest wind driven waves at moderate fetches (nearing full development) fit a consistent line on a 
plot of nondimens ional roughness versus inverse wave age. 

We also have identified candidate variables that are responsible for the large scatter on such plot for 
other field cases. We also have suggestions for why laboratory and field data do not seem to be 
relatable by a single fit on such a plot without reasonable scatter. Each category of experimental data 
are much better explained by individual fits of positive slope. The difference, we believe, is due to 
1 ) the limited ability of the laboratory facilities to simulate atmospheric turbulence, secondary flows 
involving the whole planetary boundary layer and the turbulence in the water; and 2) variability of 
currents. In other words, the state of the fluid on both sides of the interface is important for surface 
wave development. 

We therefore encourage experimenters to carefully include at least some measure of these additional 
parameters in future work. A list of the parameters needed for the air and water boundary layers and 
the wave field and, the parameters defining the context of the measurements are found in Table I. 

The need to pursue simpler conclusive experiments for practical applications is abundantly clear in 
view of the remote sensing needs, for engineering purposes and for coupled atmosphere-ocean 
models. 
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Table I. Boundary Layer, Interfacial and Context Quantities to be Measured or Estimated for a 
Compete Wind-Wave Field Experiment 


Near Surface. Local Parameters : Symbols : 

Turbulent fluctuations, mean wind velocity u’, v’, w’, U, 

and wind stress (at two or more heights, preferably) (t, u*) 

Air temperature and humidity at known height, T a , q, 

sea surface temperature, stratification parameter T s , z/L 

Directional wavenumber spectrum over all wavelengths . <J>(k) 

Slope or significant slope of waves ak, § 


Wave breaking, incidence, type 

Surface currents and their turbulent fluctuations u w , u w \ v w \ w vv * 


Context Parameters: 

Fetch F 

Duration T 

Water depth or oceanic mixed layer depth D, Z m 

Height of atmospheric PBL Z; 

Secondary flow structures 
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